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> ■ Abstract 

■ It has recently been suggested that the quark-gluon plasma formed in heavy-ion collisions behaves 

O ; 

. as a nearly ideal fluid. That behavior may be understood if the quark and antiquark mean-free- 



paths are very small in the system, leading to a "sticky molasses" description of the plasma, as 
advocated by the Stony Brook group. This behavior may be traced to the fact that there are 
p, • relatively low-energy qq resonance states in the plasma leading to very large scattering lengths for 

r~| , the quarks. These resonances have been found in lattice simulation of QCD using the maximum 

^ entropy method (MEM). We have used a chiral quark model, which provides a simple representation 

}^ • of effects due to instanton dynamics, to study the resonances obtained using the MEM scheme. In 

the present work we use our model to study the optical potential of a quark in the quark-gluon 
plasma and calculate the quark mean-free-path. Our results represent a specific example of the 
dynamics of the plasma as described by the Stony Brook group. 
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I. INTRODUCTION 



The description of the quark-gluon plasma in terms of hydrodynamics has been advo- 
cated by the Stony Brook group [ 1-3 ] . That description appears to be in accord with the 
experimental data. In such a description the motion of the quarks is characterized by an 
extremely short mean-free-path. The origin of that behavior is thought to be due to the 
relatively low-energy resonances in the qq system leading to very large scattering lengths. 
These resonances have been found in lattice studies of QCD which make use of the maximum 
entropy method (MEM) [4-9]. Similar resonances are found in the scalar, pseudoscalar, vec- 
tor and axial- vector qq channels [10]. Recently, an extensive exploration of charmonium 
studies in the confined and deconfined regions using lattice methods has been reported in 
Ref. [11]. In that work results are given for the dependence of the resonance excitation on 
the total momentum of the qq pair. We have studied that dependence for light quark systems 
in Ref. [12] and have found similar behavior to that reported in Ref. [11]. (We will make use 
of the results presented in Ref. [12] in the present work in which we calculate the imaginary 
part of the optical potential and the mean-free-path for a quark in the quark-gluon plasma.) 
We use a chiral model with rather large momentum cutoff. That model is meant to provide 
an approximate description of the instanton dynamics advocated by the Stony Brook group 
[1-3]. Earlier work using our model may be found in Refs. [13-15]. 

In our studies of meson spectra at T = and at T < Tc we have made use of the Nambu- 
Jona-Lasinio (NJL) model. The Lagrangian of the generalized NJL model we have used in 
our studies is 



C = q{ij - m')q + ^ Y^iq^qf + [qil.X'qf] (1.1) 

i=0 

- ^E[(^A^7Mg)' + (^A^757.9)^] 

i=0 

+ ^{det[g(l + As)?] + det[g(l - Ag)?]} + C,onf 

Here, m° is a current quark mass matrix, m° — diag[m^,m%w!l). The Aj are the Gell- 
Mann (flavor) matrices and A° = •\/2/31, with 1 being the unit matrix. The fourth term is 
the 't Hooft interaction and Cconf represents the model of confinement used in our studies 
of meson properties. 
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In the study of hadronic current correlators [13-15] it is important to use a model which 
respects chiral symmetry when mP — 0. Therefore, we make use of the Lagrangian of Eq. 
(1.1), while neglecting the 't Hooft interaction and Cconf- In order to make contact with the 
results of lattice simulations we use the model with the number of flavors = 1. Therefore, 
the A* matrices in Eq. (1.1) may be replaced by unity. We then use 

C = m - m')q + + [qiimf] (1-2) 

in order to calculate the hadronic current correlation functions. Thus, there are essentially 
three parameters to consider, G5, Gy and a Gaussian cutoff parameter a, which restricts the 
momentum integrals through a factor exp[— k /ot^]- As suggested by the Stony Brook group, 
we consider the NJL model and the associated chiral Lagrangian of Eq. (1.2) as providing a 
simplified representation of the instanton dynamics important for the problems considered in 
this work. Since the results obtained for the hadronic current correlation functions are similar 
in the scalar, pseudoscalar, vector and axial- vector channels, we carry out our calculations 
for the scalar qq states and multiple our results for the optical potential by 4. The parameters 
G and a were fixed in our earlier studies [12]. We take G =1.0 GeV^^ and a = 4.4 GeV. 
These values provide good fits [12] to the hadronic current correlation functions found in 
the lattice studies [10]. In order to calculate the optical potential for a quark we consider 
the quark moving in an antiquark distribution characterized by a temperature-dependent 
occupation factor n{pi) which depends upon the chemical potential The energy of a 
quark is given by E{~p ) = ["j?^ + m^]^/^. We follow the work of Shuryak [1], for example, 
and put 171 = 1 GeV. In Shuryak's work this mass is not the current quark mass, but is 
called the "chiral mass". (We would prefer to call the 1 GcV mass, the "thermal mass", 
however, the terminology used is not important for this work.) The quark thermal mass is 
given in Ref. [16], with = 4/3, as 

m^ = l^^CHT^ + 4), (1-3) 

8 TT^ 

for the case of a finite chemical potential. The thermal gluon mass is 

^l = l9"T\CA + lNf). (1.4) 

(The relation between thermal masses in QED and QCD is given on p. 146 of Ref. [16].) 
In studies of baryon matter, the chemical potentials used are often about 300 MeV or 
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less. However, once we introduce a thermal mass of about 1 GeV, we need to determine the 
chemical potential for the quarks. In this work we will consider a chemical potential of about 
1 GeV, although the calculations are easily made for other values. Once we put m = 1 GeV, 
the chemical potential is the only parameter which is varied in our study. The organization 
of our work is as follows. In Section II we discuss the calculation of the imaginary part of 
the quark optical potential and the quark mean-free-path. Section III contains some further 
discussion and conclusions. (Appendixes A and B provide a description of the calculation of 
correlators in our model. That description is readily taken over to obtain a representation 
of the qq scattering matrix. The calculation of the vacuum polarization function J(P°, P) 
which appears in Section II is described in the appendices.) 

II. CALCULATION OF THE QUARK OPTICAL POTENTIAL 

It is useful to contrast the calculation of the quark optical potential with the calculation of 
the nucleon optical potential that is to be used in the Dirac equation. The latter calculation 
is discussed in detail in Ref. [17]. In that calculation of the nucleon- nucleus potential one 
calculates the T-matrix for nucleon-nucleon scattering using the one-boson-exchange (OBE) 
model. In that case the mesons of the OBE model undergo t-channel and u-channel exchange 
between the nucleons. The result is that the imaginary part of the optical potential has a 
magnitude of about 10 MeV [17]. That in turn leads to a mean-free-path of about 10 fm 
for a 500 MeV nucleon. Many years ago, the relatively large value for the nucleon mean- 
free-path lead to the characterization of the optical model for nucleon-nucleus scattering 
as the "cloudy crystal ball" model. When we study quark propagation in the quark-gluon 
plasma we may consider a similar calculation of the optical potential at finite temperature. 
In the case of the quark-antiquark interaction the T matrix is dominated by s-channel 
resonances of the type found in the MEM studies. As we will see, the interaction in this 
case is quite strong, leading to a small mean-free-path. The resulting model is called the 
"sticky molasses" model [1-3] as opposed to the "cloudy crystal ball" model used to describe 
nucleon-nucleus scattering. We now consider the potential seen by a quark of momentum 
P2 and average over the quark spin S2- (We will consider quarks of a single fiavor, since 
that was done in the MEM studies that we have used to fix the parameters of our model.) 
In Ref. [17] the relativistic optical potential was denoted as Xl( p*, s) and, for this work, we 
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consider 

= ^^u{p^,S2)^P^,S2)u{p^,S2). (2.1) 

It is useful to introduce [14] 



2 

S2 



Here the factor of TV = 4 takes into account the sum of the interactions in the scalar, 
pseudoscalar, vector and axial-vector channels which are taken to be equal for the purposes 
of this work. The approximate equality of the interactions in these channels may be seen in 
Ref. [10]. (Note that values of Ui^p) are given in Ref. [17] for the case of nucleon- nucleus 
scattering.) 

If pi is the momentum of the antiquark in the medium, we may introduce the four- vector 

P^={pi+P2Y. (2.3) 

Now 

P2 = - P2 (2.4) 

= (E(^)+E(^))2-(^2^^2^2piP2COS^). (2.5) 

Here, we take along the z axis. We define 



HPl,P2) = 



G 



i-Gj{p,,r2)v ^^-^^ 

where J(pi , P2) is the qq vacuum polarization function defined in Appendix B. (We remark 

that we may also use the notation i(P^,p2) for the quantity defined in Eq. (2.6).) 
It is also useful to introduce the occupation factor n(pi): 

with /? = l/T and E{^i) = [^^ ^ ^2|i/2_ ^^^^^^ 

Uip2, S2)uip2, S2) = (^^) , (2.8) 

S2 

Y,^{p,,s,)v{p„s,) = , (2.9) 



5 



and note that 



Tr 



1^2 + rn\ f j) 



2m 



m 



2m 



{E1E2 -pl-pl)-m^ 



m^ 



(2.10) 



Thus, 



1 f dpl_ 
"27 (27r)= 



m 



Im 



G 



i-Gj(ri,P2) 



X 



E1E2 —pi-p2—rri^ 



(2.11) 



and 



m 



dpi 



m^ 



m 



2 E(p^)7 (27r)3E(pt 



(2.12) 



X 



E1E2 -pi-p2-m^ 



m^ 



n{pl). 



Here Ei — E(pl), E2 — E{p^) and we have made use of Eqs. (2.2) and (2.6). Values 
of i(P^,p2) are shown in Fig. 1 for values of \p^\ ranging from 0.01 GeV to 0.31 GeV. In 
Fig. 2 we show the values of n{pl) for the three values of // considered here and in Fig. 3 
we present values of lmU{p^) for those values of In Fig. 4 we show the values for the 
mean-free-path 

1 



A 



m ImU{p2) 



(2.13) 



III. DISCUSSION 

Information is available concerning the baryon chemical potential. That chemical poten- 
tial is parameterized in Ref. [18] as 

1270 MeV 

I^B = 75=7> (3-1) 

and varies strongly with ^/Snn, which is in GeV units in Eq. (3.1). For \/SViv — 200 GeV, 
we have /x^ = 26.7 MeV. 

Of particular significance for our results is the choice of the chemical potential for the 
quarks. We have taken // ~ 1 GeV. In the case of the quarks, the choice of // ~ 1 GeV leads 
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P' (GeV') 



FIG. 1: Values of t{P^,p2) are shown for various values of the quark momentum \p2\. Starting with 
the uppermost eurve, the \p2\ values in GeV units are 0.01, 0.03, 0.05, 0.07, 0.09, 0.11, 0.13, 0.15, 
0.17, 0.19, 0.21, 0.23, 0.25, 0.27, 0.29 and 0.31. (For large we have t{P^,P2) ^ {l/7rP^)G.) 
Here = {pi +P2)^, where pi is the antiquark momentum. 




FIG. 2: Values of n{pi) are shown for jj, = 1.1 GeV (dotted curve), /j, = 1.3 GeV (dashed curve) 
and II = 1.5 GeV (solid curve). Here T = 1.5 with = 270 MeV. 
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FIG. 3: The imaginary part of the quark optical potential is shown for = 1.1 GeV (dotted curve), 

fi = 1.3 GeV (dashed curve) and /j, = 1.5 GeV (solid curve). (We recall that the nucleon-nucleus 
imaginary optical potential is about 0.01 GeV in magnitude [17].) 




FIG. 4: Values of X{p2) are shown for /x = 1.1 GeV (dotted curve), /j, = 1.3 GeV (dashed curve) 
and /X = 1.5 GeV (solid curve). 
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to small mean-free-paths consistent with the suggestion of Shuryak that the resonances seen 
in the MEM analysis of the lattice results are responsible for the small mean-free-paths of 
the "sticky molasses" model. 

We remark that in nuclear matter the baryon density is 0.17 fm^'^ or 0.51 quarks/fm^ if 
we consider the nucleon to be composed of three quarks. For the values of n{pl) shown in 
Fig. 2 for the case /i — 1.3 GeV, we may calculate the density of antiquarks to be 5.91 fm~^, 
so that the density of quarks and antiquarks is about 12 fm~^ in our model. We remark 
that the energy density at RHIC for ^/S]^ = 200 GeV is 4.1 GeV/fm^ [18], which is about 
26 times the energy density of nuclear matter, which is approximately 0.16 GeV/fm^. 

In this work we have attempted to provide a quantitative analysis of the suggestion [1] 
that the large qq resonant scattering cross sections are responsible for the small quark mean- 
free-paths, with the associated relevance of the hydrodynamic description of the system that 
is created in high-energy nucleus-nucleus collisions. We have considered the interaction in 
the scalar, pseudoscalar, vector and axial-vector channels of the qq system. It is possible 
that there are important resonances of qq character, as well as the qq resonances considered 
here. Such qq states are depicted in Fig. 8a of Ref . [1] . It would be of interest to see if such 
states are found in lattice studies using the MEM scheme. 

APPENDIX A 

For ease of reference, we present a discussion of our calculation of hadronic current 
correlators taken from Ref. [ 15 ]. The procedure we adopt is based upon the real-time finite- 
temperature formalism, in which the imaginary part of the polarization function may be 
calculated. Then, the real part of the function is obtained using a dispersion relation. The 
result we need for this work has been already given in the work of Kobes and Semenoff 
[19]. (In Ref. [19] the quark momentum is k'^ and the antiquark momentum is k'^ — P^. 
We will adopt that notation in this section for ease of reference to the results presented in 
Ref. [ 19 ] .) With reference to Eq. (5.4) of Ref. [ 19 ] , we write the imaginary part of the scalar 
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polarization function as 

x{[l - ni{k) - n2(A;)]5(F° - Ei{k) - E^ik)) 
-[ni{k) - n2(A;)]5(P° + Er{k) - E^^k)) 
-[n2{k) - ni{k)]S{P' - Ei{k) + E^ik)) 
-[1 - n^{k) - n2(A;)]5(P° + E^{k) + E^ik))} . 

Here, Ei{k) = [P + ml{T)Yl'^ . Relative to Eq. (5.4) of Ref.[19], we have changed the sign, 
removed a factor of g'^ and have included a statistical factor of Nc- In addition, we have 
included a Gaussian regulator, exp[— fc^/a^]. The value a = 0.605 GeV was used in our 
applications of the NJL model in the calculation of meson properties at T = 0. We also 
note that 

^i(^) = eW) + l ' ^^2) 

and 

^^(^) ^ ePE,l) + 1 • (A3) 

For the calculation of the imaginary part of the polarization function, we may put k"^ — 
m\{T) and {k — Py — m|(T), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (Al) the factor 13s arises from a trace involving Dirac matrices, such 
that 

= -Tr[(^ + mi)(^-f + m2)] (A4) 
= -2(mi + m2)^ (A5) 

— * 

where mi and m2 depend upon temperature. In the frame where P = 0, and in the case 
mi = m2, we have Ps = 2Po^(l — Aw? / Pq). For the scalar case, with mi = m2, we find 

Im r, = ^ (l - i^)'" e-''"-'[l - 2nm ] , (A6) 



where 



^' = ^-m'(r). (A7) 
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For pseudoscalar mesons, we replace (5s by 

(3p = -Tr[z75(|^ + mi)275(^-f + ^2)] (A8) 
= - 2(mi - m2f , (A9) 

which for mi = 1712 is (3p = 2Pq in the frame where P — 0. We find, for the tt mesons, 

Im MP', r) = ^ (1 - e-' 1 _ 2n,(k) ] , (AlO) 

where k"^ = Pq/4 — m^{T), as above. Thus, we see that the phase space factor has an 
exponent of 1/2 corresponding to a s-wave amphtude. For the scalars, the exponent of the 
phase-space factor is 3/2, as seen in Eq. (A6). 
For a study of vector mesons we consider 

Pi = + - f + ^2) ] , (All) 

and calculate 

^A[P^ -ml-ml + 4mim2 ] , (A12) 

which, in the equal-mass case, is equal to 4Po + Sm'^{T), when P — 0. This result is 
needed when we calculate the correlator of vector currents. Note that, for the elevated 
temperatures considered in this work, m„(T) = rridiT) is quite small, so that 4Pq + 8m^(T) 
can be approximated by 4Pq , when we consider the vector current correlation functions. In 
that case, we have 

Im Jv{P\ T) ~ ^Im Jp(p2, T) . (A13) 
At this point it is useful to define functions that do not contain that Gaussian regulator: 



,2 ^,_iVcPoV. 4m^(T)V/^ 



Im Jp{P\ T) = ^ ( 1 - ] [ 1 - 2ni(fc) ] , (A14) 



and 

Im MP\ T) \^ {l - l!^) ''V - 2m(^) ] , (A15) 

For the functions defined in Eq. (A14) and (A15) we need to use a twice-subtracted dispersion 

relation to obtain Re Jp(P^,r), or Re Jy(P^,T). For example, 

p2 . 

Re Jp(p2, T) = Re Jp(0, P) + ^[ Re Jp(Po , T) - Re Jp(0, T) ] (A16) 

"0 

P\P'-Pi) Im Jp(.,T) 



JAm'^iT) S{P^ - S){P^ - S) ' 
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where can be quite large, since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce Jp{P'^,T) and Jv{P'^,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re Jp(P^,T), or Re Jv'(P^,T), by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
Jp{P^T) and MP^T). 

In order to make use of Eq. (A16), we need to specify Jp(0) and Jp{Pq). We found it 
useful to take P^ = -1.0 GeV^ and to put Jp(0) = Jp(0) and Jp{P^) = Jp{Pq)- The 
quantities Jv(0) and Jv{Pq) are determined in an analogous function. This procedure in 
which we fix the behavior of a function such as ReJy(P^) or ReJv'(P^) is quite analogous 
to the procedure used in Rcf. [ 20 ] . In that work we made use of dispersion relations to 
construct a continuous vector-isovector current correlation function which had the correct 
perturbative behavior for large P^ — > — oo and also described the low-energy resonance 
present in the correlator due to the excitation of the p meson. In Ref. [ 20 ] the NJL model 
was shown to provide a quite satisfactory description of the low-energy resonant behavior 
of the vector-isovector correlation function. 

We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

^C(P^ T) = j d^xe^^-^ « T{j{x)m) » , (A17) 

where the double bracket is a reminder that we arc considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp^i{x) = 
q{x)ij5X'q{x), where, in the case of the tt mesons, i — 1,2 and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i{x) — q{x)X^q{x), where i = for the flavor-singlet current 
and i = 8 for the flavor-octet current. 

In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 
of the basic vacuum polarization function of the NJL model, Jp{P^,T). Thus, 

where Gp{T) is the coupling constant appropriate for our study of tt mesons. We have 
found Gp{T) = 13.49 GeV""^ by fitting the pion mass in a calculation made at T = 0, with 
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TTT'u — nT'd — 0.364 GeV. The result given in Eq. (A18) is only expected to be useful for small 
P^, since the Gaussian regulator strongly modifies the large behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P^. 



Cp{P',T) 

p2 



Jp{P\T) 

p2 



l-Gp{T)Jp{P\T)- 



(A19) 



(As usual, we put P = 0.) This form has two important features. At large Pq, 
lm.Cp{PQ,T)/ Pq is a constant, since Im Jp(Po,T) is proportional to P|. Further, the 
denominator of Eq. (A19) goes to 1 for large Pq. On the other hand, at small Pq , the 
denominator is capable of describing resonant enhancement of the correlation function. As 
we have seen, the results obtained when Eq. (A19) is used appear quite satisfactory. ( We 
may again refer to Ref. [20], in which a similar approximation is described.) 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
J/i,i(3^) — Q.{x)liJ,\(l{x) with i=l, 2 and 3. In this case we define 



and 



(A20) 



(A21) 



taking into account the fact that the current jii,i{x) is conserved. We may then use the fact 
that 



Jv(p^^) = ^^,.J^(p^^) 



(A22) 



and 



Im Jy(P^T) 



P^ + 2ml{T) 
Stt 



^ImJp(P^^). 



1 - ^^^^ e-^^^/"^ [ 1 - 2n,{k) ] (A23) 



p2 



(See Eq. (A7) for the specification oi k = \k\.) We then have 

1 



Cv(P\T) ^ MP\T) 



l-Gv{T)MP\T) 



where we have introduced 



ImJy(P^r) 



P2 + 2ml{T) 



Stt 



~ -ImJp(P^^). 



4ml{T) 



p2 



1/2 



2ni{k) 



(A24) 

(A25) 

(A26) 
(A27) 
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In the literature, uj is used instead of Pq [ 4-6 ] . We may define the spectral functions 

cry (a;, T) = ^ Im Cy (c^, T) , (A28) 

TT 

and 

ap{uj, T)^-lm Cp{uj, T) , (A29) 

TT 

Since different conventions are used in the literature [4-6], we may use the notation 
CTp(a;, T) and av{uj, T) for the spectral functions given there. We have the following relations: 

ap{uj,T)^ap{uj,T), (A30) 

and 

av{uj,T) 3 

^ = -l<^v{^, T) , (A31) 

where the factor 3/4 arises because, in Refs. [4-6], there is a division by 4, while we have 
divided by 3, as in Eq. (A22). 



APPENDIX B 

— * 

Here we extend the work of Appendix A to consider case of finite three-momentum, P. 
We consider the calculation of ImJp(P°, P, T). The momenta P° and P are the values 
external to the loop diagram. Internal to the diagram, we have a quark of momentum 
k'^ -\- P^ jl leaving the left-hand vertex and an antiquark of momentum — P'^/2 entering 
the left-hand vertex. It is useful to define 



Ei{k) = k + P/2 



p2 \ 1/2 

k'^ + — + kPcose] 



(Bl) 
(B2) 



and 



E2{k) = k-P/2 



= {k^ + — -kPcos9 



1/2 



(B3) 
(B4) 



Here k 



and P = \P\ 
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We have 



Here, 



and 



ImJp(P°,P,r) 



27r 



2Ei{k)2E2{k) 



x{[l - mik) - n2(A;)]5(P° - E,{k) - E^ik)) 
-Mk) - n2{kMP' + Ei{k) - E^ik)) 
-[n2{k) - ni{k)]S{P'' - E,{k) + E^^k)) 

-[1 - m{k) - n2(A;)]5(P° + E^{k) + E2{k))} . 



ni{k) = 



ef^Mk) + 1 ' 



n2{k) 



(B5) 



(B6) 



(B7) 



In Eq. (B5), the second and third terms cancel and the fourth term does not contribute. It 
is useful to rewrite 5(P° — Ei{k) — E2{k)) using 

2 



where 



5[/(cos^)] 



df 



-5{cos6 — x) , 



I d cos d I X 



(B8) 



We find 



and obtain 





-- cos'^e 






_ AP^ie + PV4) - Po^ 






4^2p2 


df 


2 


Ei{k)-E2{k) 


9 cos 6' 


Ei{k)E2{k) 



lmJp{P',P,T) = ^N,Ppe{P'){27rf J 



k'^dk 



I 



2Ei(k)E2(k) 



[1 - ni{k) - n2ik)] 



(27r)3 
df(cose) 



dcos 9 
x5{cos9 — x)d{cos6) . 



(B9) 



(BIO) 



(Bll) 
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We note there is a singularity when Ei{k) = E2{k). That occurs when cos 6* = or 6 = 7^/2. 
For our calculations we eliminate the point with 6 = tt/2 when evaluating the angular 
integral over (i(cos6')5(cos6' — x) in the last expression. We obtain 



ImJp(pO,P^r)=iV,/3pe(P°) 



An' 



^max 



(27r)3 

[l-n^{k)-n2{k)] 



X 



kP\Ei{k) - E2{k)\ 



where x is obtained from Eq. (B9), 



X 



pO 

kP 



4 4 



1/2 



(B12) 



(B13) 



For the calculations reported in this work we have P^ = E(pl) + E(p^) and P = ih+P2, 
where p2 is the quark momentum and pi is the antiquark momentum. Thus, we may also 
use the notation J{pi,P2) as we have done in the main text. 
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